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Abstract. The extra-solar planet HD 209458b has been found to have an extended atmosphere of escaping atomic 
hydrogen (Vidal-Madjar et al. 2003), suggesting that "hot Jupiters" closer to their parent stars could evaporate. 
Here we estimate the atmospheric escape (so called evaporation rate) from hot Jupiters and their corresponding 
life time against evaporation. The calculated evaporation rate of HD 209458b is in excellent agreement with the 
Hi Lyman-Q observations. We find that the tidal forces and high temperatures in the upper atmosphere must 
be taken into account to obtain reliable estimate of the atmospheric escape. Because of the tidal forces, we show 
that there is a new escape mechanism at intermediate temperatures at which the exobase reaches the Roche lobe. 
From an energy balance, we can estimate plausible values for the planetary exospheric temperatures, and thus 
obtain typical life times of planets as a function of their mass and orbital distance. 
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1. Introduction 

Among the more than one hundred extra-solar planets 
known, over 15% are closer than 0.1 AU from the central 
star. Except for a possible transiting planet (OGLE-TR- 
56b, Konacki et al. 2003; Torres et al. 2003), there are no 
detection of planets closer than about 0.04 AU. Following 
the discovery that the planet HD 209458b shows a sur- 
prisingly large escape of atomic hydrogen (Vidal-Madjar 
et al. 2003; confirmed by Vidal-Madjar et al. 2004), the es- 
cape flux from hot Jupiters needs to be evaluated. Already 
in 1995, together with the release of the first hot Jupiter 
discovery (Mayor & Queloz 1995), Burrows & Lunine 
(1995) highlighted the evaporation issue. But Guillot et 
al. (1996) concluded that the mass loss was not signifi- 
cant. However their estimates did not consider two factors 
which need to be addressed: the infiuence of the strong 
tidal forces from the neighboring parent stars, and the 
high upper atmosphere temperature (Hebrard et al. 2003). 

2. The upper atmosphere model 

Radiative equilibrium is commonly used to calculate the 
temperature of the upper atmosphere (see Schneider et al. 
1998). But this is not appropriate because it does not ap- 
ply to the low density upper atmosphere. As an example, 
in the Solar system, the temperature of planetary upper 
atmospheres (thermosphere, exosphere) is much higher 



Send offprint requests to: A. Lecavelier des Etangs, e-mail: 
lecaveliSiap . f r 



than Tcff, the effective temperature of the lower atmo- 
sphere. For Earth and Jupiter, Tog are ~250K and 150 K, 
while thermospheric temperatures are ~1000K for both 
planets (Chamberlain & Hunten 1987). Although the ob- 
served high temperatures in the giant planets are not yet 
explained, extreme and far ultraviolet fluxes. Solar wind 
and perhaps gravity waves may contribute to the heating 
in uncertain amounts {e.g., Hunten & Dessler 1977). 

The atmospheric escape flux strongly depends on the 
temperature profile of the atmosphere. As the tempera- 
ture of the upper atmosphere of extra-solar planets is not 
known, in a first step, we will use it as an input parameter: 
Tup. We consider a simple atmospheric structure similar 
to that observed in the giant planets of the Solar system, 
that is a two level temperature structure: Toff in the lower 
atmosphere and Tup above. We assume that the lower at- 
mosphere (up to the thermobase) is mainly composed of 
molecular hydrogen at Tcs. Above that level, in the ther- 
mosphere and exosphere, the gas is a mixture of atomic 
and molecular hydrogen at Tup (Fig.nj. 

The quantitative results presented in this letter do 
not depend on the estimated Toff, but on the density 
at the thermobase, rifji and rifja for atomic and molec- 
ular hydrogen, respectively. In the following we will use 
three different estimates of these densities. In the case of 
Jupiter and Saturn, H2 is the main constituent at the ther- 
mobase, where the temperature transition is observed to 
occur for nni — 2 x 10^ cm^'^ and — 1 x 10^^ cm^'^ 
(Chamberlain & Hunten 1987); these values will be used 
for the Model A. By evaluating the production of atomic 
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Fig. 1. Vertical structure of the atmosphere of 
HD 259458b as a function of the temperature of the 
upper atmosphere (Model A). At the top of the thermo- 
sphere, the exobase is the critical level where the mean 
free path is equal to the distance to the Roche lobe. For 
high temperatures, the exobase reaches the Roche lobe; 
this leads to a geometrical blow-off. 



hydrogen with realistic condition for the lower atmosphere 
of HD 209458b, Liang et al. (2003) have shown that the 
Hi mixing ratio could be about 10% at the thermobase. 
We thus use the corresponding nni = 1 x 10^° cm~^ and 
riHa — 9x 10^° cm~^ for the intermediate Model B. Finally, 
at high thermospheric temperatures, it is likely that most 
of the molecular hydrogen is dissociated into atomic hy- 
drogen (Coustenis et al. 1998); we use nni = 2x 10^^ cm~'^ 
and riHa = cm~'^ for the corresponding Model C. 



3. The atmospheric escape 

The atmospheric escape critically depends on atmospheric 
temperatures. Below a critical temperature, Tc, the escape 
flux can be calculated using the Jeans escape estimate 
which refers to the escape of particles whose velocity is 
in the tail of the Boltzmann distribution and that have 
enough energy to escape the planets gravity. For tem- 
peratures above Tc, the kinetic energy of the atoms or 
molecules is sufficient to overcome gravitational forces and 
they can stream, or blow-off. For Jeans escape the flux is 
calculated at a critical level corresponding to the exobase 
(Hunten et al. 1989). This critical level is the level above 
which particles can freely escape without collisions. The 
exobase is usually defined by the place above which the 
mean free path (l/nQ) is larger than H, the scale height 
of the atmosphere, where n is the volume density. This 
definition results from the approximate calculation of the 
integrated collisional cross-section (Q) from the exobase 
to infinity: J^^i,^^^ n{r)Qdr « Uc^oHQ = 1, where 
the density at the exobase. Here this idea must be gener- 
alized to take into account the particular geometry of hot 
Jupiters; we define the exobase by the place above which 
the mean free path is larger than the distance to the Roche 



lobe of the planet: f^°f^'^''°^'^ n(r)Qdr — 1. The Roche 

^ J exobase ^ ' ^ 

lobe is the last equipotential around the planet beyond 
which the equipotentials are open to infinity or to encom- 
pass the star. Thus the exobase is the level above which 
atoms and molecules can definitively escape the planet. 
Note that in the case of an isolated planet, because the 
Roche lobe is at infinite distance, the two definitions are 
identical. Then, the total escape flux from the planet is 



given by M — An 
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/i is the mass of the escaping elements (Hi and H2), rexo 
is the radius of the exobase, and vt is the thermal veloc- 
ity at temperature T (Chamberlain & Hunten 1987). A 
is defined by A = — xm/^^j where x is the gravitational 
potential. We calculated the density profile in the atmo- 
sphere using the barometric law: Vn = nVA. Recapture 
of escaped particles is not possible because of the high ra- 
diation pressure pushing away hydrogen atoms at hundred 
of kilometer per seconds (Vidal-Madjar et al. 2003). 

We can see from Fig-^that between the classical Jeans 
escape at low temperatures and the dynamical blow-off of 
the atmosphere at high temperatures, the tidal forces lead 
to a new escape mechanism at intermediate temperatures. 
Indeed, for HD 209458b (d=0.047 AU, Mp=0.69Mjup, and 
i?p=1.35Rjup) and Tup above ~5000K, the exobase comes 
close to the Roche lobe, and the kinetic energy is similar 
to the potential energy needed to reach this limit. In that 
case, the escape mechanism is a geometrical blow-off which 
is due to the filling up of the Roche lobe by the thermo- 
sphere pulled up by the tidal forces. Indeed, the dynamical 
blow-off is obtained when the escape velocity Wosc is of the 
order of (or smaller than) the thermal velocity vt-, that is 
for small value of Acxo ('^— ^osc/^t)- With the characteris- 
tics of HD 209458b, the dynamical blow-off takes place for 
temperatures above 20 000 K. However, at temperatures 
between 5 000 K and 20000 K, although Aoxo ^1, we have 
Acxo— ARocho ^1, and most of the gas at the exobase can 
freely reach the Roche lobe. This geometrical blow-off is 
due to the spatial proximity of the Roche lobe through 
which the gas can escape the planet. 

4. Tidal forces 

In the case of hot Jupiters, the planet's gravity is sub- 
stantially modified by stellar tidal forces. The common 
assumption has been to neglect tidal forces by considering 
isolated planets far from their star. But tidal forces have a 
significant influence on the density distribution in the up- 
per atmosphere of hot Jupiters. In order to calculate the 
potential energy x, we include the difference between the 
stellar gravitation and the centrifugal effect in the orbiting 
planet reference frame, which results in the tidal forces. 
The equipotentials are modified, from quasi-spherical in 
the lower atmosphere to asymmetric elongated shapes at 
the level of the Roche lobe. As a result, Xi hence the verti- 
cal density distribution, strongly depend on the location in 
both longitude and latitude on the planet {6, cj)). We calcu- 
late the escape rate as the sum of M{d, cj)), the escape flux 
as a function of longitude and latitude. Even with a uni- 
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Fig. 2. Shape of the exobase and the corresponding es- 
cape rate (Model A, Tup=11100K). Light grey is for the 
large escape rates (10^° gs~^ str~^) and black is for the 
small escape rates (<5xl0^gs~^ str~^). The star is on 
the X-axis toward positive coordinates. The escape rate is 
the largest toward the star and in the opposite direction. 



form temperature, the escape flux per unit area is larger 
toward the star and in the opposite direction (Fig. 

Finally, using the above, we estimate the total escape 
flux of atomic and molecular hydrogen as a function of the 
upper atmosphere temperature for HD 209458b (Fig. 
For temperature above 8 000 K (Model A), 7 000 K (Model 
B) and 6 000 K (Model C), the Hi escape flux is larger 
than the minimum flux of 10^*^ gs~^ needed to explain 
the occultation depth of 15% as observed in Lyman a 
(Vidal-Madjar et al. 2003). From the escape rate, we can 
derive the corresponding life time needed to evaporate the 
total mass of the planet ti = AIp/M. However when the 
planet mass decreases, the evaporation rate increases. This 
results in a shorter life time given by i2 = / dM/M, which 
is typically shorter than ti by a factor of 5 to 10 (Fig.^D). 

5. The temperature of the upper atmosphere 

Although the mechanism responsible for the heating of the 
upper atmospheres in the Solar system is not fully iden- 
tified, we obtain a plausible estimate of Tup from a com- 
parison of heating and cooling mechanisms. A lower limit 
of the heating can be estimated from the energy flux of 
both the stellar extreme ultraviolet (EUV) and Lyman a 
photons. Using the EUV and a monte-carlo simulation of 
the multiple scattering of Lyman a photons in the ther- 
mosphere, we estimate the minimum energy input to be 
G = 2.7 X (d/l AU)^^ ergcm^^ s^^, where d is the orbital 
distance to the star. The cooling is due to a combina- 
tion of the heat conduction toward the cooler lower at- 
mosphere, coUisional exitation of the Hi electronic levels, 
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Fig. 3. Escape flux and life time of HD209458b as a func- 
tion of Tup. The results from various models show the 
uncertainty due to the lack of constraint on the densities 
at the thermobase. 

(a) shows the escape flux of Hi (long dashed) and H2 
(short dashed) in the Model A, and the total escape flux 
(solid lines) for models A, B and C. For the Model A, the 
total escape flux from an isolated planet (without tidal 
forces) is plotted with the dot-dashed line. The tidal forces 
increase the evaporation rate by one to two orders of mag- 
nitude. 

(b) shows the planet life time. A first estimate is obtained 
from the instantaneous escape flux for Model A (ii, dashed 
line). Taking into account the planet's mass decrease with 
time, we find shorter life times (^2, solid lines). 
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Fig. 4. Upper atmosphere temperature estimated from 
the energy balance as a function of the orbital dis- 
tance of a 0.69 Mjup planet. The input energy is 2.7 x 
(d/lAU)"^ ergcm^^ s^^. Below 0.04 AU, the temperature 
decreases because the large escape rate cools down the 
thermosphere. 
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Fig. 5. Contour plot of the planet life time (logj^g *2/year) 
as a function of the mass and orbital distance (Model A) . 
The temperature is obtained from the energy balance and 
is varying with time. We obtain similar numerical esti- 
mates and contour plots for Model B and C. The squares 
show the positions of detected planets. OGLE-TR-56b is 
at the upper left Conner. HD 209458b is shown by a cir- 
cle. HD 49674b, HD 46375b, HD 76700b, and HD 83443b 
(triangles) have life time of about 10^ years. These plan- 
ets must have lost a large fraction of their hydrogen and 
might be remnants of former hot Jupiters. 

collisional ionization (photo- ionization is negligible), and 
cooling by escaping atoms and molecules carrying off their 
kinetic energy. We derive the ionization fraction and the 
coUisonal cooling following Spitzer (1978) and Osterbrock 
(1989) for atomic and molecular hydrogen clouds. With 
the HD 209458b characteristics, the ionization fraction is 
barely above 0.01. We calculated the heat conduction fol- 
lowing the derivation given by Hunten & Dessler (1977). 

If we apply this calculation to Jupiter with UV heat- 
ing only, we find a temperature rise at the thermobase as 
low as ~15 K (Strobel & Smith 1973); clearly other heat- 
ing mechanisms are at work in Jupiter (Hunten & Dessler 
1977). Considering the UV heating only, the energy bal- 
ance provides a lower limit to the thermospheric tempera- 
ture (Fig. HI). For HD 209458b, we evaluate the lower hmit 
of Tup to be about 1 1 100 K, 9 800 K and 7 900 K for Models 
A, B and C, respectively. Cooling is dominated by atmo- 
spheric escape and collisonal excitation of H i. With these 
temperatures, the present evaporation rates of atomic and 
molecular hydrogen from HD 209458b are estimated to 
be Mm ^ 5.7 x 10i°gs-i and M^^ ~ 1.7 x lO^Ogs-^ 
in Model A, Aim ^ 1-3 x 10" gs'^ and ^ 3.5 x 

lO^gs"! in Model B, and Mm ^ 5.2 x 10" gs^^ in 
Model C, in agreement with the observational lower limit 
of ~10^°gs-i (Vidal-Madjar et al. 2003). We obtain a 
life time (i2) of 10^" to 10^^ years. During 5 x 10^ years, 
HD 209458b may have lost 1% to 7% of its initial mass. 

With models of heat conduction within escaping at- 
mospheres, Watson et al. (1981) showed that the escape 
rate can be limited by heat exchange. Using their result 
applied to HD 209458b, we found that the escape rate is 



heat-limited to ~10^^ gs ^, in agreement with the results 
obtained by Trilhng (1999) and Lammer et al. (2003). 
Although these results neglect the tidal forces, our val- 
ues are below this limit and therefore do not need a full 
calculation of the heat exchange within the thermosphere. 

At a given temperature, the escape rate is roughly pro- 
portional to the thermobase densities and the energy in- 
put. However, using the temperature calculated from the 
energy balance, the escape rate is less sensitive to the as- 
sumed input parameters. It is roughly proportional to the 
thermobase densities to the power of 0.3, and proportional 
to the energy input to the power of about 0.4. 

6. Life time of hot Jupiters 

Finally, we can also estimate the life time of a given 
planet as a function of its mass and orbital distance 
(Fig. . We conclude that planets with orbital distance 
lower than 0.03-0.04 AU (corresponding to orbital periods 
shorter than 2-3 days) have short life time unless they 
are significantly heavier than Jupiter. This may explain 
why only few planets have been detected with periods be- 
low 3 days. Low-mass hot Jupiters have also short life 
times, meaning that their nature must evolve with time. 
These planets must loose a large fraction of their hydro- 
gen. This process can lead to planets with an hydrogen- 
poor atmosphere ("hot Neptunes"), or even with no more 
atmosphere at all. The emergence of planets modified by 
evaporation (and possibly the emergence of the inside core 
of former and evaporated hot Jupiters) may constitute a 
new class of planets (see also Trilling et al. 1998). If they 
exist, these planets could be called the "Chthonian" plan- 
ets in reference to the Greek deities who come from hot 
infernal underground (Hebrard et al. 2003). 
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